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1  Introduction 

Direct-detection  (DD)  optical  fiber  frequency  division  multiple  access  (FDMA)  is  a  simple  and  practical  alter¬ 
native  to  optical  heterodyne  FDMA.  Previous  works  on  estimating  the  performance  of  DD  optical  FDMA  networks 
have  usually  focused  only  on  the  (linear)  crosstalk  degradation,  typically  relying  on  simplifying  approximations  such 
as  the  dominance  of  the  adjacent  channel  interference  [1],  [4].  Other  sources  of  performance  degradation  such  as 
signal  loss  due  to  optical  (predetection)  and  electrical  (postdetection)  filtering,  intersymbol  interference  (ISI)  due  to 
the  optical  filtering,  and  channel  beats  (or  nonlinear  crosstalk)  are  usually  ignored. 

This  paper  presents  a  more  precise  analysis  that  takes  into  account  the  effect  of  optical  and  electrical  filtering, 
ISI,  and  linear  crosstalk  [3].  The  model  used  here  is  valid  for  arbitrary  optical  filter  transfer  functions  and  received 
pulse  shapes.  We  consider  in  this  paper  only  On-Off-Keyed  (OOK)  modulation;  extending  the  analysis  to  Frequency- 
Shift-Keyed  (FSK)  modulation  is  straightforward,  but  is  not  discussed  here. 

2  System  Model 

The  system  is  assumed  to  consist  of  M  transmitter-receiver  pairs  connected  by  a  star  coupler.  We  assume  for 
convenience  that  M  is  even.  The  received  signal  r(t)  at  each  receiver  is  given  by 

M/2—1  oo 

r(t)  =  V 2  X  X  diks{t  -  kT)cos(2wfit  +  <j>i)  (1) 

i—~M/2k  =  ~oo 

where  dik  {0, 1}  is  the  data  bit  transmitted  by  channel  i  during  the  fcth  bit  period  kT  <  t  <  (k  +  1)T,  and 
<f>i  are  the  carrier  frequency  and  phase  of  channel  i  respectively,  and  s(t)  is  the  received  pulse  shape,  normalized  to 
unit  energy.  Note  that  all  M  channels  are  assumed  to  be  bit-aligned.  This  results  in  a  worst  case  linear  crosstalk  as 
shown  in  [2].  Laser  phase  noise  and  frequency  chirp  are  ignored.  We  choose  channel  0  (the  center  channel)  to  be  the 
signal,  with  the  other  M  —  1  channels  being  crosstalk.  Strictly  speaking,  this  is  with  some  loss  of  generality  since 
the  filter  frequency  response  is  a  weak  function  of  the  chosen  channel.  However,  for  most  cases,  this  effect  is  small 
and  can  be  neglected. 

The  front-end  optical  filter  (see  Fig.  1)  is  assumed  throughout  this  paper  to  be  linear  time  invariant  (LTI)  with 
an  impulse  response  h(t)  and  with  a  system  bandwidth  of  P  Hz.  For  optical  filters  with  an  aperiodic  frequency 
response,  P  equals  the  filter  tuning  range;  for  periodic  optical  filters  P  equals  the  frequency  response  period,  or  “free 
spectral  range”.  The  channel  spacing  Af  equals  P/M.  Shot  noise  is  assumed  negligible  in  comparison  to  the  receiver 
thermal  noise,  as  is  usually  the  case  in  practical  direct  detection  receivers.  The  photodetector  thus  acts  as  an  ideal 
square-law  detector  with  output  current 

t  2 

bi(n  =  Rr<  (  [  b(r)r(t  -  r)dr )  (2) 

Wo  /  LP 

where  RPd  is  the  photodetector  responsivity  (henceforth  set  to  1)  .  The  subscript  LP  denotes  removal  of  the  double 
frequency  terms. 

The  total  receiver  current  is  i(t )  — -  i<t{t)  +  4,(1),  where  the  noise  current  in(t)  is  assumed  to  be  a  zero  mean 
AWGN  process  with  double  sided  spectral  density  No/2.  The  current  i(t)  is  subsequently  low  pass  filtered  (modeled 
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by  an  integrate  and  dump  filter),  sampled,  and  compared  to  a  threshold.  The  sampled  output  for  the  Oth  bit  period 
is  given  by 

Z(d,M)  =  f  id(t)dt  + 

Jo 

The  variable  £  is  the  received  energy,  and  is  a  function  of  the  M  data  bits  d  =  {di,  i  =  —M/2, ...  ,  M/2  —  1} 
sent  by  the  M  channels  during  [0,T].  The  noise  variable  Af  is  a  zero  mean  Gaussian  random  variable  with  variance 
<t2  =  NqT/2.  Assuming  negligible  ISI,  we  may  approximate  the  received  energy  during  each  bit  period  by  that 
collected  with  an  infinite-time  integrator 


/ 


in{t)dt  —  £(d,  M)  +  Af. 


(3) 


/OO  pOO 

id(t)dt=  /  \H(f)\2\R(f)\2df. 

'OO  J  —  OO 


(4) 


This  form  simplifies  the  analysis  considerably.  Later,  we  shall  correct  for  this  approximation. 


The  energy  £(d,  M)  is  a  sum  of  three  components:  a  signal  component  <S(do,  M ),  a  linear  crosstalk  component 

C(d~ ,  M ),  and  a  nonlinear  crosstalk,  or  channel  beats,  component  B( d,  M ).  Here,  by  abuse  of  notation,  d~  =  d  —  do- 
The  three  components  of  £  can  be  written  as 


/OO 

|F(/)|2|5(/-/0)|2d/ 

-OO 


(5) 


l^i 

The  sampled  output  is  then 


M/2-1 

C(d~  ,M)  =  di 

i=z-M/2 
0 

\H(f)\2S(f  -  +  \H(f)\2 S(f  +  +  ft)e^-^df. 


Z{ d,  M)  =  <S(d0,  M)  +  C(d~  ,M)  +  B(d,  M)  +  Af 


/OO 

\H{f)\2\S{f  ~  h)\2df 

-OO 


(6) 

(7) 

(8) 


i.e.,  the  signal  is  corrupted  by  two  non-Gaussian  noises  in  addition  to  the  Gaussian  noise  source.  The  expressions  in 
(5)-(7)  will  enable  us  to  determine  the  crosstalk  degradation  and  signal  energy  reduction  (due  to  the  optical  filter). 
This  is  done  in  the  following  Section.  Later,  we  shall  determine  the  additional  losses  caused  by  the  electrical  filtering 
(due  to  the  finite  time  integration)  and  the  ISI. 


3  Analysis 


Define  coo 

£i{di,M)^di  /  |if(/)|2|5(/-/*)|2d/  -,di  g  {0,1} 

J  — OO 


(9) 


as  the  received  energy  from  channel  i  when  it  sends  di.  Since  the  filter  is  centered  on  f0,  the  carrier  frequency  of  the 
signal,  we  may  write  (9)  as  f00 

£i(di,  M)  =  di  J  \H(6)\2\S(S-±)\2d6  (10) 

where  S  ==  f/P  is  normalized  frequency.  In  particular,  the  signal  energy  equals 

5(d0,  M)  =  £0(d0,  M)  =  d0  I  |H'(«5)|2|5(<5)|2dl5.  (11) 

J  —  OO 


Further  simplification  depends  on  whether  |.ff(<5)|2  is  an  periodic  or  aperiodic  function.  In  the  former  case  we  may 
convert  to  the  time  domain  to  obtain 

OO 

£i(di,M)  =  di  J2  %]%)e-i2,r?i/^ 

q—  —  oo 


(12) 


where  /i[q]  are  the  Fourier  coefficients  of  |if(<5)|2  and  s(n)  are  integer  samples  of  the  Fourier  transforms  of  |S(<5)|2. 
In  particular,  the  signal  energy  when  a  ONE  is  sent,  5(1,  M ),  equals 

OO 

5(1,  M)  =  E  %]s(g).  (13) 

q~  —  oo 


Note  that  we  assume  an  ideal  extinction  ratio,  which  gives  £(0,  M)  =  0.1  The  received  energy  is  thus  maximized 
when  all  M  channels  send  ONE: 

M/2-1  oo 

£max(M)  =  E  E  %]S(g)e-2'2XS,'/M 

fc  — M/2  ?=- oo 

OO 

=  Af  E  h[<lM]s(qM).  (14) 

q—  —  o o 

From  (13)  and  (14),  the  maximum  crosstalk  energy  equals 

OO  OO 

Cmam(M)  =  £mam(M)-S(l,M)  =  M  E  h[qM]s{qM)~  E  %]%)•  (15) 

q~  —  oo  qzz—oo 


Similar  expressions  are  obtained  for  optical  filters  with  an  aperiodic  transfer  function.  A  summary  of  these  expressions 
are  given  in  Table  1. 

From  (12)  above  we  may  also  obtain  the  received  energy  probability  mass  function  (pmf),  which  is  given  by  the 
convolution  of  the  M  independent  random  variables 


£i(di,  M) 


fO  ;  di  =  0  prob  =1/2 

(  £i(l,M)  ;  d{  =  1  prob  =  1/2, 


where  —M/2  <  i  <  M/2  —  1,  giving 


Pe{£) 


M/2-1 


i  =  — M/2 


±s(e)  +  ±6(e-£i{i,M))  , 


(16) 


(17) 


where  0  denotes  convolution  and  S(-)  is  the  Dirac  delta  function.  The  crosstalk  pmf,  denoted  as  C(d~ ,  M),  is 
obtained  by  removing  the  i  =0  term  in  (17). 


3.1  Filtering  Penalty,  Crosstalk  Penalty  and  BER  Computation 

Ideally,  the  (normalized)  received  signal  energy  when  a  ONE  is  sent  is  1.  Filtering  the  signal  energy  reduces  the  signal 
energy  to  5(1,  M).  Therefore,  the  signal  filtering  penalty  equals 


Xfu  =  10  log 


10 


5(1, M)‘ 


(18) 


Similarly,  the  crosstalk  penalty  is  the  ratio  of  the  signal  distance  with  worst-case  crosstalk  to  that  without  crosstalk. 
The  former  is  equal  to  5(1,  M)  —  Cmax(M),  while  the  latter  is  5(1,  M).  The  crosstalk  penalty  is  thus 


X, 


10  log 


5(1,  Af) 


10 


5(1,  Af )  —  Cmax(M) ' 


(19) 


To  determine  the  bit  error  rate  (BER)  degradation,  we  first  determine  the  expected  energy  output  conditional 
on  the  data  sent  by  the  crosstalk  channels 

E(Z\dv  =  0,d-)  =  C(d“,M) 

E(Z\do  =  1,  d_)  =  5(1,  M)  -f  C(d“ ,  Af) 


1  Extending  the  analysis  to  include  non-zero  extinction  ratios  is  simple,  but  not  considered  here. 


where  d~  is  defined  in  Section  2.  The  noise  variance  is  equal  to  <r2  =  NqT/2  for  both  ONE  and  ZERO.  Thus, 
conditional  on  knowledge  of  the  crosstalk  pmf  C(d~,  M),  our  problem  becomes  the  known  problem  of  detecting 
binary  signals  in  AWGN,  with  the  threshold  7  equal  to  (<5(1,  M)  +  Cmax(M))/ 2.  This  gives  a  probability  of  error 


M 


5(1,  M)  +  C(d~,  Af)  —  7 


(20) 


where  the  averaging  is  performed  over  the  crosstalk  pmf  C(d  ,  M)  and  Q(x)  =  -i-e  X*l2dx. 

4  Performance  of  FP  filters 

4.1  Frequency  and  Time  Responses 

An  FP  filter  with  mirror  (power)  reflectivity  R  has  a  (field)  impulse  response  h(t)  given  by 

OO 

h(t)  =  (1  -  R)  RiHt  -  i/P) 

and  intensity  frequency  response 


i=o 


W)i2  =  rrl  E 


The  Fourier  coefficients  of  |ff(/)|2  are  thus 

h[n] 


R 


f?l,: 


00  <  n  <  00. 


1  +  2? 

Alternatively,  the  squared  frequency  response  is  given  by  the  Airy  function 

=  1  +  (2F/n)2  sin2(^) 


(21) 

(22) 

(23) 

(24) 


where  the  Finesse  F  =  tt^/R/{1  —  R)  and  the  Free  Spectral  Range  (or  period)  P  —  c/2nL,  with  n  the  refractive 
index  inside  the  cavity.  The  FP  filter’s  full  width  at  half  maximum  (FWHM)  B  equals  B  =  P/F. 

Fig.  2  shows  the  pulse  response  of  an  FP  filter  to  a  pulse  with  F  —  100  and  1/BT  =  0.5  at  resonance  ( 6  =  0) 
and  slightly  off  resonance  (6  =  0.03).  Notice  the  “relaxation  resonance”  in  the  off-resonance  pulse  response. 

We  find  it  convenient  here  to  define  the  ratios 


a  =  1/BT  =  data  rate  /  filter  FWHM 
(3  =  1/TAf  =  data  rate  /  channel  spacing. 

Note  that  a  is  filter  dependent,  while  f3  is  filter  independent.  For  the  FP  filter  a  =  F/TP  while  /3  =  M/TP ,  giving 
(3 /a  =  M/F  (=  filter  FWHM  to  channel  spacing  ratio). 

4.2  Penalty  Calculations 

For  a  received  NRZ  pulse  shape,  s(t)  =  1/ \/T  for  t  £  [0,  T]  and  zero  elsewhere.  It  follows  that 


|5(<5)|2  =  TP 


[wTP6)\ 


kTPS 


J 


whose  transform  equals 


»(r) 


1  -  r<TP 

0  r  >  TP. 


(25) 

(26) 


Substituting  (23)  and  (26)  in  (12)  and  (14), and  assuming  F  ^  1  and  [3  <  1  we  obtain 

an 


5(1,  M) 


1  —  ^rjCosech2(7r/22’)  «  1  —  a/x 


(27) 


£max  (  M  ) 


JV 


irM 

~2F 


zot)y{*M/2F) 


-  cosech2  (  kM  /2F ) 


(28) 


where  we  used  the  approximations  coth(s)  «  l/x  +  x/3  and  cosech  (a)  m  1/x  —  ®/6,  which  are  valid  for  small  x.  It 
follows  from  above  that  the  worst-case  crosstalk  is  given  by 


Cm**{M)  =  £max{M)-S{l,M) 


7t 

12 


(29) 


The  above  formulas  agree  precisely  with  our  previous  results  in  [2],  where  we  assumed  a  =  0  (i.e.,  zero-bandwidth  or 
tone-like  channels).  In  that  case  the  signal  energy  and  the  crosstalk  are  given  by  1  and  ^  (^)  respectively.  Fig.  3 
shows  the  energy  pmf  when  an  FP  filter  is  used  with  M/F  =  0.8.  The  two  distinct  and  identical  masses  represent 
the  possible  energy  levels  received  when  the  signal  sent  ZERO  or  ONE.  This  figure  also  shows  the  effect  of  increasing 
a,  which  is  to  be  interpreted  here  as  increasing  the  data  rate.  The  left  (zero)  mass  shifts  to  the  right  (increasing 
crosstalk  energy),  while  the  right  (one)  mass  shifts  to  the  left  (increasing  the  signal  energy  loss).  Both  effects  result 
in  a  closing  of  the  eye-opening  (equivalently,  a  reduction  in  signal  distance). 


By  substituting  (27)  in  (18)  we  obtain  the  signal  filtering  penalty 

1 


Xfu  =  101og10 


1 


(30) 


which  is  shown  in  Fig.  4  (bottom  solid  line)  as  a  function  of  a.  However,  as  mentioned  earlier  in  Section  2,  not  all 
the  filtered  signal  energy  is  collected  by  the  integrate  and  dump  filter  since  the  integration  is  only  over  a  finite  period 
of  length  T  and  not  over  all  time.  This  loss  is  referred  to  here  as  the  finite  time  integration  penalty.  For  general  low 
pass  filters,  it  is  the  electronic  filtering  penalty.  In  the  case  of  the  FP  filter,  optimizing  the  integration  period  (of 
duration  T )  of  a  received  pulse  yields  the  signal  energy  1  -  which  is  slightly  less  than  1  -  a/rr.  Thus,  the 

finite  time  integration  penalty  is 


Xfin  =  10  log 


10 


,1  -  (In  2  +  1/2)2/  '  (31) 

The  third  signal  related  penalty  is  the  ISI  penalty  which  is  a  result  of  detecting  energy  from  neighboring  pulses  in 
addition  to  the  desired  pulse.  By  considering  only  nearest  neighbor  pulses,  we  find  that  the  ISI  in  turn  reduces  the 
signal  distance  from  1  —  (In  2  +  1/2)2  to  1  —  (2 In  2)2,  yielding  an  additional  penalty  equal  to 


Xisi  —  10log10 


/l  — (ln2+l/2)g\ 
V  1  —  (2  In  2)2  j 


(32) 


Finally,  we  have  the  crosstalk  penalty.  Here  we  use  the  signal  energy  after  both  optical  and  electrical  filtering  as  well 
as  ISI  as  the  baseline  (or  numerator  of  log  argument),  giving  a  crosstalk  penalty  equal  to 


Xcr  —  10  log10 


(2  In  2)2 


1  —  (21n2)2  _  £~(M/_F)2(1 


(33) 


Fig.  4  is  a  plot  of  the  cumulative  penalty  as  a  function  of  a.  The  first  three  penalties  above  are  signal-related,  and 
thus  are  zero  for  a  =  0.  The  crosstalk  penalty,  however,  is  nonzero  even  for  a  =  0,  due  to  the  M/F  factor. 


4.3  BER  degradation 

The  BER  degradation  with  increasing  a  is  obtained  by  (numerically)  computing  the  crosstalk  pdf  C(d~,M),  and 
using  it  as  shown  in  (20).  An  example  is  shown  in  Fig.  5  for  M/F  =  0.4.  Here  the  signal  to  noise  ratio  (SNR)  is 

defined  as  SNR  =  1/cr,  where  <r  =  N0T/2  is  the  noise  variance.  Thus  the  SNR  is  referred  to  the  input  of  the  optical 
filter. 


5  Performance  of  Multistage  MZ  filters 


5.1  Physical  Structure  and  Power  Transfer  Function 


An  ideal  single-stage  Mach-Zehnder  (MZ)  filter  consists  of  two  ideal  3-dB  couplers  connected  by  two  waveguides  of 
length  difference  A L  (see  Fig.  6(a)).  The  power  transfer  function  of  the  two  output  ports  referenced  to  the  same 


input  port  are  given  by  ,  , 

|tfi(/)f=cos2(^)  \H2{f)\2=^(~) 


(34) 


where  P  =  c/nAL  is  the  transfer  function  period. 

To  construct  a  multistage  MZ  filter,  single-stage  MZ  filters  are  cascaded  in  a  manner  such  that  the  length  difference 
of  stage  i  is  double  that  of  stage  i  —  1,  with  the  first  stage  having  a  length  difference  A L  as  in  Fig.  6(b)  (see  also  [5], 
[6]).  The  overall  transfer  function  of  a  If -stage  MZ  filter  is  given  by 


K 

\H{S)\Z  =  JJcos2(7r2i-15)  = 

i=l 


sm2(nMS) 
M 2  sin2(?r<5) 


M  =  2k 


(35) 


with  Fourier  coefficients 


AT— |n| 
0 


|n|  <  M 
otherwise. 


(36) 


Fig.  7  shows  the  transfer  functions  of  three  stages  comprising  a  3-stage  MZ  filter  and  Fig.  8  shows  one-period  of 
the  resulting  transfer  function.  By  placing  all  but  the  signal  channel  in  the  2K~l  nulls  of  one  period  of  the  overall 
transfer  function  and  the  signal  channel  in  the  passband,  the  crosstalk  may  be  strongly  suppressed.  However,  the 
non-zero  channel  bandwidths  as  well  as  any  filter  imperfections  naturally  prevent  a  total  crosstalk  suppression.  The 
filter  FWHM  B  is  narrowed  as  the  number  of  stages  increases,  since  for  large  K  we  may  show  that 


B 


2.8  P 


The  envelope  of  the  transfer  function  is  closely  approximated  by  an  Airy  function  (Eq.  (24)) 


(37) 


\H(6)\2 


1 

1  +  M 2  sin2(ir<5)  ’ 


where  the  Finesse  F  =  tM/2.  The  accuracy  of  this  approximation  can  be  seen  from  Fig.  9.  Note  that  this  Airy 
function  has  a  smaller  FWHM  than  that  of  the  MZ  filter  it  approximates.  Thus  an  FP  filter  with  the  above  Finesse 
will  produce  both  more  crosstalk  and  more  signal  loss  than  the  MZ  filter.  As  a  result,  the  performance  of  this  FP 
filter  is  a  lower  bound  on  the  MZ  filter  performance.  For  example,  an  ideal  MZ  filter  with  K  —  7  stages  will  perform 
no  worse  than  an  FP  filter  of  finesse  F  &  200. 


5.2  Temporal  Responses 


By  construction,  the  ith  stage  of  a  /f-stage  MZ  filter  has  two  paths  of  length  difference  A  Li  —  2t~lAL,  where  A  L 
is  the  length  difference  of  the  first,  or  “coarsest”,  stage.  In  effect  then,  an  input  impulse  “sees”  exactly  2K  different 
paths  of  lengths  0,  AA, . . . ,  (2K  —  1)AL.  Thus  a  A'-stage  MZ  filter  is  a  finite  impulse  response  (FIR)  filter  with  an 
impulse  response  given  by  a  finite  train  of  impulses  of  spacing  At  —  nAL/c  —  1/P  and  of  equal  height  2~K  =  1/M. 

In  contrast  to  our  analysis  of  the  FP  filter  model,  we  shall  find  it  more  convenient  to  express  our  results  in  terms 
of/3  instead  of  a.  However,  these  two  quantities  are  dependent  variables  here  since,  by  using  (37), 


B_ _ 1_ 

A  J  BT 
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Fig.  10  shows  the  pulse  response  of  a  7-stage  MZ  filter  to  two  (unit  intensity)  pulses  with  (3  —  0.25  and  2.  As  can 
be  seen,  for  minimal  pulse  distortion  /3  has  to  be  smaller  than  1,  allowing  the  pulse  response  to  be  concentrated  in 
a  T-second  long  period,  and  to  approach  its  ideal  value  of  1  during  that  period. 


5.3  Penalty  Calculations 

Substituting  (36)  and  (26)  in  (12),  we  obtain 
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The  last  equation  is  exact  and  is  independent  of  f3.  It  follows  directly  from  (39)  and  (40)  that 

5(1,  M)  ta  1  —  3/3  for  large  M 

and  a  i 


3  M2  -  1 

Cf nax(M)  =  -( — — )  ss  3/3  for  large  M. 


The  filtering  penalty  is  thus 


Xfu  —  10  Iog10 


i-3/z)- 


Finding  the  finite  time  integration  penalty  is  done  as  with  the  FP  filter;  the  integration  period  position  is  chosen  to 
maximize  the  received  energy.  Some  simple  algebra  reveals  that  the  maximum  received  energy  during  a  T-second  is 
1  —  6/5/12,  thus  giving  a  finite  time  integration  penalty  of 
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Similarly,  the  ISI  penalty  is  found  following  the  same  procedures  used  with  the  FP  filter  to  be 
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From  above,  the  crosstalk  penalty  becomes  independent  of  M  for  large  M.  Again,  by  measuring  the  penalty  fiom 
the  filtering+finite  time  integration+  ISI  signal  energy  baseline,  we  get 

=  10Io«‘»  {^m)  ■  <«> 

Fig.  11  is  a  plot  of  the  resulting  cumulative  penalty. 

5.4  BER  Degradation 

Computing  the  BER  degradation  follows  steps  identical  to  those  used  in  the  FP  filter  example.  The  difference  is  in 
the  values  of  5(1,  M )  and  C(d~ ,  M).  Fig.  12  shows  the  BER  degradation  with  3  {or  K  —  A  stages. 


6  Multistage  MZ  Filter  Imperfections 

In  this  Section  we  discuss  two  crucial  imperfections  unique  to  multistage  MZ  filters.  These  imperfections  may 
cause  the  performance  obtained  with  these  filters  to  be  seriously  degraded,  and  are  therefore  worth  considering. 

0.1  Imperfect  Coupling  Ratios 

An  ideal  multistage  MZ  filter  is  constructed  with  ideal  3  dB  couplers.  If  the  coupling  ratio  deviates  from  50:50,  then 
the  filter  response  changes  significantly.  Let  |ei;|,  |t2i|  <1  be  the  deviations  of  the  coupling  ratios  of  the  first  and 
second  3  dB  coupler  in  the  fth  MZ  stage.  Using  a  first  order  series  approximation,  we  find  that  the  crosstalk  penalty 
with  5  =  0  is  given  by 

10logl  + 


1  2  5Zt=l  lli 


where  yu  —  e$-  4-  e\{  and  72  i  =  (ei;  —  f2i)2 3 4 5 6 ■  The  filtering  plus  crosstalk  penalty,  on  the  other  hand,  depends  only  on 
the  71;  and  equals 

10  logic  f.  * K - )•  (48) 
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This  penalty  is  shown  for  several  values  of  Ky  in  Fig.  13,  where  we  have  assumed  for  simplicity  that  yu  =  7. 

0.2  Stage  Misalignment 

Another  possible  imperfection  in  MZ  filters  is  an  error  in  the  misalignment  of  the  stages  such  that  the  resonances  do 
not  coincide  as  in  Fig.  7.  Using  a  second-order  series  approximation  and  after  considerable  algebra,  we  find  that  the 
filtering  plus  crosstalk  penalty  may  be  expressed  as 

101°8l°  G-2*2£f=1€?)  (49) 

where  e;  is  the  normalized  deviation  of  stage  i  from  its  ideal  length,  the  normalization  factor  being  the  signal 
wavelength  \o  —  c/  fo-  Fig.  14  shows  the  increase  in  crosstalk  energy  due  to  stage  asynchronism  for  several  values  of 
K  and  the  sensitivity  increase  with  the  number  of  stages. 

7  Summary  and  Conclusions 

We  have  presented  an  analysis  of  the  performance  of  direct  detection  optical  fiber  FDMA  networks  in  the  presence 
of  linear  crosstalk,  optical  and  electrical  filtering,  and  signal  ISI.  We  used  the  developed  model  and  theory  to 
analyze  the  performance  of  FP  filters  and  multistage  MZ  filters  as  channel  demultiplexers.  Our  results  show  that 
the  degradations  due  to  signal  filtering  and  linear  crosstalk  are  roughly  of  equal  magnitude  for  practical  system 
operating  points,  with  finite  time  integration  and  ISI  effects  being  of  less  significance.  However,  filter  imperfections 
such  as  filter  loss  (not  discussed  here),  nonideal  coupling  ratios,  and  stage  misalignment  may  significantly  reduce 
the  obtainable  performance.  The  additional  effect  of  channel  beats  is  complicated  by  the  need  to  assume  channel 
bit  -asynchronism  and  will  be  considered  in  future  research. 
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Figure  11:  A  plot  of  the  various  penalties  with  an  MZ 
filter  as  a  function  of  (3  for  large  K.  Notation:  F:  filter¬ 
ing.  F+FI:  filtering  +  finite  time  integration.  F+FI+ISI: 
filtering  -f  finite  time  integration  +  ISI. 


Figure  13:  Effect  of  splitting  ratio  mismatch  on  the 
crosstalk  penalty  of  a  6-stage  MZ  filter  with  OOK.  Here 
Tit  =  T2i  =  7- 


Figure  12:  BER  degradation  with  an  8-stage  MZ  filter 
and  OOK. 


Figure  14:  Increased  crosstalk  energy  due  to  stage  asyn- 
chronism  as  a  function  of  e.  Sensitivity  increases  with 
number  of  stages,  m  is  the  order  of  the  resonance  used 
for  channel  selection 
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Figure  6:  (a)  Ideal  single  stage  MZ  filter  with  length  d- 
ifferential  A L  (b)  Construction  of  a  K-stage  MZ  filter. 
Each  stage  is  marked  with  the  appropriate  length  differ¬ 
ential. 


Figure  7 :  Frequency  responses  of  each  stage  of  a  3-stage 
MZ  filter  Overall  transfer  function  of  the  3-stage  MZ  fil¬ 
ter. 


Figure  9:  Approximating  the  upper  envelope  curve  of  an 
MZ  frequency  response  with  a  frequency  response  of  an 
“equivalent”  FP  filter.  Case  shown  is  for  if  =  10  and 
F  =  1608. 
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Figure  10:  Pulse  response  of  a  7-stage  MZ  filter  with 
normalized  input  frequency  5  =  0.  Shown  in  the  figure 
are  the  response  due  to  two  pulse  of  duration  /3  <  1  and 

f3  >  1. 


Figure  8:  Overall  transfer  function  of  a  3-stage  MZ  filter. 
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Figure  1:  Receiver  structure:  Optical  filter  -  Detector  - 
Electric  low  pass  filter.  p(<  T )  is  the  time  shift  required 
to  maximize  the  received  signal  energy. 
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Figure  4:  A  plot  of  the  various  penalties  for  different  M/F 
with  an  FP  filter  as  a  function  of  a.  Notation:  F:  filter¬ 
ing.  F+FI:  filtering  +  finite  time  integration.  F+FI+ISI: 
filtering  +  finite  time  integration  +  ISI.  Other  curves  in¬ 
clude  F+FI+ISI  plus  crosstalk  at  several  values  of  M/F. 


Figure  2:  Pulse  response  of  an  FP  filter  with  F  =  100  to 
a  pulse  of  width  T  =  200/P. 
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Table  1:  Comparison  of  Aperiodic  and  Periodic  Filter 
Output  Energy  Expressions 
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Figure  5:  BER  crosstalk  induced  degradation  with  an  FP 
filter  with  M/F  =  0.4.  Also  shown  for  comparison  is  the 
BER  curve  for  one  channel  (no  crosstalk). 


Figure  3:  Energy  pmf  with  FP  filter.  P=100,  M— 80. 


